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We study the persistent current of correlated spinless electrons in a continuous one-dimensional 
ring with a single weak link. We include correlations by solving the many-body Schrodinger equation 
for several tens of electrons interacting via the short-ranged pair interaction V(x — x'). We solve 
this many-body problem by advanced configuration-interaction (CI) and diffusion Monte Carlo 
(DMC) methods which rely neither on the renormalisation group techniques nor on the Bosonisation 
technique of the Luttinger-liquid model. Our CI and DMC results show, that the persistent current 
(/) as a function of the ring length (L) exhibits for large L the power law typical of the Luttinger 
liquid, I oc L" 1 "", where the power a depends only on the electron-electron (e-e) interaction. For 
strong e-e interaction the previous theories predicted for a the formula a = (1 + 2c(rg) 1 ^ 2 — 1, where 
ctRG = [V(0) — V(2kF)]/2irhvF is the renormalisation-group result for weakly interacting electrons, 
with V(q) being the Fourier transform of V(x — x'). Our numerical data show that this theoretical 
result holds in the continuous model only if the range of V(x — x') is small (roughly d < 1/2&F, 
more precisely 4d 2 k F <C 1). For strong e-e interaction (org > 0.3) our CI data show the power law 
/ oc L~ x ~ a already for rings with only ten electrons, i.e., ten electrons are already enough to behave 
like the Luttinger liquid. The DMC data for org > 0.3 are damaged by the so-called fixed-phase 
approximation. Finally, we also treat the e-e interaction in the Hartree-Fock approximation. We 
find the exponentially decaying I(L) instead of the power law, however, the slope of log(/(L)) still 
depends solely on the parameter org as long as the range of V(x — x') approaches zero. 

PACS numbers: 73.23.-b, 73.61.Ey 



I. INTRODUCTION 

A clean one-dimensional (ID) wire biased by contacts 
with negligible backscattering is known to exhibit the 
conductance quantized as an integer multiple of e 2 /h. 
The effect can be explained in a Fermi-liquid model of 
non-interacting quasi-particles?^ In fact, a clean ID sys- 
tem is not a Fermi liquid due to the electron-electron (e-e) 
interaction. Away from the charge-density-wave instabil- 
ity, the system is a correlated Luttinger liquid with collec- 
tive bosonic elementary excitations which contrast with 
independent fermionic quasi-particles of ordinary Fermi 
liquids-. Nevertheless, the Luttinger-liquid model gives 
for a clean ID wire the same conductance (e 2 /h per spin) 
as the Fermi-liquid model 4 *. 

If a localized scatterer is introduced into the wire, 
the conductance quantization breaks down. For non- 
interacting electrons, the Landauer formula^ expresses 
the conductance per spin as (e 2 /h)\tk F \ 2 , where tk F is 
the transmission amplitude through the scatterer at the 
Fermi level. 

In the Luttinger liquid models, the infinite wire which 
contains a single structureless scatterer, exhibits the con- 
ductance varying with temperature as oc T 2a , where a 
depends only on the e-e interaction. Thus, for a > (re- 
pulsive e-e interaction) the wire is impenetrable at T — > 
regardless the strength of the scatterer. A similar power 
law exhibits at T — > the differential conductance as a 



function of the bias voltage (oc U 2a ) and conductance 
versus the wire length (oc L~ 2a ). These power laws are a 
sign of the Luttinger liquio£&£. The T 2a and U 2a laws 
were successfully measured^. 

In this work we deal with similar power laws exhib- 
ited by interacting electrons in an isolated mesoscopic ID 
ring. In particular, magnetic flux <f> piercing the opening 
of the mesoscopic ring gives rise to the persistent electron 
current circulating along the ringi. This persistent cur- 
rent is given at T = OK as / = — dEQ(4>) / d<j), where Eq is 
the energy of the many-body groundstatc. If the ring is 
clean and the single-particle dispersion law is parabolic, 
the e-e interaction does not affect the persistent current 
due to the Galilean invariance of the system^. 

However, if a single scatterer is introduced, the non- 
interacting and interacting result differ fundamentally. 
For non-interacting spinless electrons in the ID ring con- 
taining a single scatterer with transmission probability 
\tk F \ 2 <C 1, the persistent current at T — OK depends on 
the magnetic flux and ring length (L) as^ 

I = (ev F /2L) \i kF \ sin(27r</>/0o), (1) 

where 0o = h/e is the flux quantum, hp is the Fermi 
wave vector, and vf is the Fermi velocity. For a spinless 
Luttinger liquid the persistent current follows the power 
law / oc More precisely^ 

/oc J L" Q " 1 sin(27T0/0o), (2) 
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where the power a depends only on the e-e interaction, 
not on the properties of the scatterer. The formulae (Q]) 
and |(2} were derived 9 - assuming large L. 

The formula can also be obtained heuristically 9 - as 
follows. Matveev et alJ£ analyzed, how the e-e interac- 
tion renormalizes the bare transmission amplitude tk F of 
a single scatterer inside the ID wire with two contacts. 
They derived the renormalized amplitude t}~ F by using 
the renormalization-group (RG) approach suitable for a 
weakly-interacting electron gas. If the system length (L) 
is large, their result can be expressed in the form 

t kp ~(t kF /\r kF \)WL) a , a>0, (3) 

where |tVf| 2 = 1 — l^fe^ | 2 , d is the spatial range of the 
e-e interaction V(x — x'), and the power a is given (for 
spinless electrons) by the expression 

a RG = [V(0) - V{2k F )]/2nhv F , (4) 

with V(q) being the Fourier transform of V(x — x'). We 
note that the formulae © and were derived assuming 

1 « \n{l/d) « a RG -\ (5) 

where I is a properly chosen scale of the RG theoryi -. 
Moreover, a = org only for ccrc "C 1 (weak e-e inter- 
action). For strong e-e interaction (say ana — 0.5) the 
theory 3,1 1 predicts the more general result, 

a = (1 + 2a RG ) 1/2 - 1. (6) 

This result is believed 1 -! to hold for any V(x — x') with 
range d which is finite but which can in principle be quite 
large (in comparison with 1//cf). Finally, if we replace in 
equation (JTJ) the bare amplitude tk F by the renormalized 
amplitude ([3]), we recover the power law where a is 
now given by the microscopic formulae ([6]) and ([4]). 

In the Luttinger-liquid model, the physics of the low- 
energy excitations is mapped onto an effective field the- 
ory using Bosonization 9 -, where terms expected to be neg- 
ligible at low energies are omitted. Within this model, 
the asymptotic dependence ^ was obtained by using 
the analogy to the problem of quantum coherence in dis- 
sipative environment 9 -. To avoid this analogy as well as 
Bosonization, in Ref^ 2 - the persistent current was cal- 
culated by solving the ID lattice model with nearest- 
neighbor hopping and interaction. Applying numerical 
RG methods, the formula ([2]) was confirmed for long 
chains and strong scatterersJ^ 

However, insofar it has not been verified whether the 
RG formula (01 holds in a microscopic model which does 
not rely on the RG approach. We present such micro- 
scopic many-body model in this work. Dealing with a 
continuous model, we can vary the range of the e-e inter- 
action in order to test the robustness of the formula (g]) 
against various shapes of V{x — x'). This point was not 
addressed in the lattice-model-based studies as the range 
of the e-e interaction was fixed to the nearest-neighbor- 
site interaction. 



Further, derivation ;] 9 ' 10 of the formulae ([2]), ([3]) and 
(HJ rely on the large number of particles limit. Here we 
directly address the interesting question what is the min- 
imum number of electrons which exhibits the onset of the 
Luttinger-liquid dependence / cx As we will illus- 

trate, such behavior can be identified from system sizes 
of the order of ten electron. 

We study the persistent current of correlated electrons 
in a continuous ID ring containing the strongly-reflecting 
scatterer, because strong backscattering is know n 9 ' 10 ' 12 
to reduce the system size necessary to achieve the L~ Q_1 
asymptotics. Using advanced configuration-interaction 
(CI) and diffusion Monte Carlo (DMC) methods, we solve 
the continuous many-body Schrodinger equation for sev- 
eral tens of electrons interacting via the e-e interaction 

V(x-x')=V exp(-\x-x'\/d). (7) 

Interaction {7} emulates screening (say by the metallic 
gates) and allows us to compare our results with the re- 
sults of the correlated models 9 - ' 10 ' 12 which also assume 
the e-e interaction of finite range. Our CI and DMC cal- 
culations do not rely on the RG techniques and solve the 
continuous model, unlike the RG studies focused largely 
on the lattice model o 12 ' 13 ' 14 ' 15 . Our numerical data show 
that the formulae (O and Q hold only if the range of 
the e-e interaction is small (d < 1/2^). For strong e-e 
interaction (ana si 0.3) our CI data show the power law 
/ cx L -1 ~ Q already for rings with only ten electrons. In 
other words, ten electrons are already sufficient to show 
the Luttinger-liquid behavior. 

It is known that the fermion sign problem causes 
an exponential inefficiency of the DMC method 1 ^ un- 
less it is circumvented by the so-called fixed-node or 
the fixed-phase approximation for inherently complex 
wave functions^. The accuracy of our DMC results for 
a RG <; 0.3 is therefore limited by the quality of the phase 
in the Hartree-Fock trial wave function which is employed 
in the fixed-phase approximation. Since Hartree-Fock is 
the simplest possible trial wave function and does not 
capture the Luttinger-liquid correlation effects it is not 
surprising that the fixed-phase bias becomes pronounced. 
We analyze these findings later in detail. (Our DMC 
should not be confused with the path-integral Monte- 
Carlo method s 18 ' 19 ' 20 ' 21 , used to study the tunnelling 
conductance of the Luttinger liquid in the bosoniscd 
model.) 

Finally, we treat the e-e interaction in the self- 
consistent Hartree-Fock calculation. We find for large 
L the exponentially decaying I(L) instead of the power 
law. However, the slope of log(/(L)) still depends solely 
on the parameter ana given by the RG formula (|4]), as 
long as the range of V(x — x') approaches zero. 

In section II. A we start with the single-particle ap- 
proach to the ID ring. In section II. B we define our 
interacting many-body model. In section II. C we outline 
how we solve the many-body model in the Hartree-Fock 
approximation. In sections II. D and II. E we describe 
how we obtain the fully-correlated many-body solution 
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by means of the DMC method and CI method. Our re- 
sults are discussed in section III and a summary is given 
in section IV. Technical details are in the appendices. 



II. THEORY 

A. Single-particle model 

We consider the circular ID ring threaded by magnetic 
flux <f> = BS = AL, where S is the area of the ring, B 
is the magnetic field (constant and perpendicular to the 
ring area), and A is the magnitude of the resulting vector 
potential (circulating along the ring circumference). In 
this section we discuss the single-particle states in such 
ring. In general, the single-electron wave functions ipn{x) 
in the ID ring obey the Schodinger equation 



fld_ + 2tt^_ r; , 

2m V i dx L <fio 



ipn(x) = e n i/}„(x) (8) 



with the cyclic boundary condition 

lp n (x + L) = 1p n (x) 



(9) 



where m is the electron effective mass, x is the electron 
coordinate along the ring, and U (x) is an external single- 
particle potential. 

We introduce the wave functions <p n (x) by substitution 



ip„(x) = ip n (x) exp ( —i^—x 
L 4> 



(10) 



If we set (|T0|) into flSJ) and ©, we obtain the equation 

ip n (x) = e n (fi n (x) (11) 



2m dx 2 



with the boundary condition 



tp n (x + L) = exp i27r— if n (x) . (12) 

Equations (| 1 1 1) and (fT2")) can be solved for an arbitrary 
potential U(x) numerically. Consider first the ring region 
x G (-L/2, L/2) as a straight-line segment of an infinite 
ID wire. Inside the segment the potential is U(x), outside 
we keep it zero. Therefore, the wave function outside is 



ip k (x) = ae lkx + be~ lkx , x < -L/2 , 

ip k (x) = ce lkx + de- ikx , x>L/2. 
The amplitudes a and b are related to c and d by 



(13) 



(14) 



(15) 



where Tq is the transfer matrix^ 



T n = 



tk 



(16) 



with tk and being the transmission and reflection am- 
plitudes of the electron impinging the region (—L/2, L/2) 
from the left. At the boundaries we express tp k (—L/2) 
and (pk(L/2) by using equations (fTB")) and (fTl)) . To come 
back to the ring threaded by magnetic flux, we relate 
ipk(—L/2) and ipk(L/2) through the boundary condition 
(|12p . Combining the obtained relation with the equations 
(fT5)l and (fTfi| we obtain the equation 



T 



^ = exp(z2-7r— ) ' 

4>o 



where 



T = 



ifei 



rk_ ~ikL J_ - 



ikL 



(17) 



(18) 



Thus exp(i27r0/0o) is the eigenvalue of the matrix T. 
The product of the eigenvalues of this matrix is given by 
its determinant which is unity. The second eigenvalue is 
thus exp(— i2n(j)/(f)o). Their sum is equal to the matrix 
traced , which gives the equation for the spectrum^, 



(19) 



H)= Rc 


exp(— ikL) 







The numerical solution of equation (fT9")) has to be com- 
bined with numerical computation of the transmission 
amplitude tk (the algorithm for computation of tk and 
is described in the Appendix A). The solution of equa- 
tion (fT5|) gives us the dependence k n ((f>) and eventually 
the single-particle eigenenergy e n (<j>) = h 2 k 2 i ((j))/2m. 

For each k n (<f>) we can also calculate the wave function 
ip n (x). We proceed as follows. By means of equations 
(fT7|) and (fT5|) we express the amplitude a in the form 



1 



tk n 



J(2ir<P/4, + k„L) 



b, 



(20) 



where the amplitude b can be obtained by normalizing 
the wave function. Then we express from equation Q13p 
the boundary conditions tp n (—L/2) and dip n (—L/2)/dx. 
Finally, we combine these boundary conditions with nu- 
merical solution of equation (jTTJ) in the discrete for 



,23 



<Pn (Xj+l) = 

2 171 ITT I \ 



(fin (Xj) - (fi n (ij-i) , (21) 



where A is the step, Xj = jA, and j — 0, ±1, ±2, 
Once the wave functions ip n (x) are known, the wave func- 
tions ip n (x) can be obtained by means of the relation (HI 
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B. Interacting many-body model 

We now consider the ring with N interacting ID elec- 
trons. This system is described by the Hamiltonian 



H 



N 

E 

3=1 L 
N 



2m \idx ' 



2tt 
L wit 



(22) 



where Xj is the coordinate of the j-th electron, j8(x) is 
the potential of the scatterer, and V(xj — Xi) is the e-e 
interaction ([7]). The eigenfunction ^{x\,X2, ■ ■ ■ , xn) and 
eigenenergy E obey the Schodinger equation 

= E^! (23) 

with the cyclic boundary condition 

$>(xi, . . . ,Xi + L, . . .,x N ) = . . . ,Xi, . . - ,x N ) (24) 

for i = 1,2, ...,N, 

If the system is in the groundstate with eigenfunc- 
tion $>o(xi,X2, ■ • ■ ,xn) and eigenenergy Eq, the persis- 
tent current can be expressed as 



where 



1 = (*nl I|*o 



a y-^ ( h <> 



mL 



dxj 



e0 

T 



(25) 



(26) 



is the TV-particle current operator. One can calculate / 
directly from the relation (|2"5|) . or one can rewrite ([23]) 
by means of the Hcllman-Feynman theorem 



dE dH . _ . 



(27) 



Using p7 | . ([26 ]) . and ([22 ]l one gets the formula^ 

I= -il E ^)- (28) 



C. Hartree-Fock approximation 

In the Hartree-Fock model, the ground-state wave 
function ^ is approximated by the Slater determinant 



*o(zi, ■ ■ -,%n) 



iPn(xi) 



ipi(x N ) 
iPn(xn) 



(29) 



The wave functions tp n (x) obey the Hartree-Fock equa- 
tion 



2 m 



9 2tt 



dx L (f> 
+ U H (x) + U F (n,x) 



j6(x) 

i> n (x) = e n tp n (x) (30) 



with the boundary condition 

lp n (x + L) = 1p n (x) 

where 



U H (x) = J2 J dx'V(x-x')\Tp n ,(x')\ 2 



(31) 



(32) 



is the Hartree potential and 
UF(n, x) = 

^J^Y, J<te'V(x-x')lP n (x'Wn'(x')Tpn>(x) (33) 



is the Fock nonlocal exchange term (expressed as an effec- 
tive potencial for further convenience). In the equations 
(|32| and (|33|) we sum over all occupied states n'. 
If we use again the substitution 



ijj n {x) = tp n (x)exp ( -i^—rx 
L 0o 



(34) 



the equations ([30" |) - ([3"3]) give the Hartree-Fock equation 



h 2 d 2 

■z — -j-x + jS(x) + U H {x) + U F (n,x) 
2m ax z 



£ n ip n (x) (35) 



with the boundary condition 

if n (x + L) = cxp i2-k- 



<Pn(x) 



(36) 



where the potentials Ujj{x) and Up(n,x) are still given 
by equations (|3"2")l and (f3"3")) , but with ip n replaced by ip n . 
From equations (|22|) - ((36| the groundstate energy E — 
(^fo\H\^fo) can be expressed as 



2 {<Pn \Uh{x) + U F {n,x)\ <p T , 



(37) 



The Hartree-Fock equation (|35[) can be solved by the 
same procedure as the single-particle equation (fTTj) as- 
suming that the potential U{x) = j5(x) + Uh(x) + 
/7i?(n, x) is known. We apply this procedure iteratively in 
order to obtain the self-consistent Hartree-Fock solution. 
In the first iteration step we solve equation (f3"5"|) for the 
non-interacting gas, i.e., for Uh(x) = and Up(n, x) = 0. 
The resulting (p n (x) is used to evaluate the Hartree and 
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Fock potentials, where the term Up(n,x) has to be eval- 
uated for each n separately. These potentials are used in 
the second iteration step to obtain new (p n (x) and new 
potentials Uh(x) and Uf(ti,x), etc., until the energies 
e„ and ground-state energy (j37|) do not change anymore. 
In the Appendix B the iteration procedure is described 
including a few nontrivial details. Setting the resulting 
ground-state energy ([57)1 into we obtain the persis- 
tent current. Of course, this Hartree-Fock calculation 
does not include the many-body correlations. To include 
the correlations we use the DMC and CI techniques de- 
scribed in the next two sections. 



D. Diffusion Monte Carlo (DMC) model 

Consider first the iV-electron ID Schodinger equation 
7i*(X) = £?tf(X) (38) 



In the DMC algorithm the equation (|41 [) is solved 
stochastically based on the action of the projection op- 
erator exp[— t(H — Et)}^- At time t = one chooses a 
proper trial wave function v&j^X). A set of walkers (sam- 
pling points in the iV-dimensional electron configuration 
space) is generated according to the distribution \ty(X). 
A small timestep r is chosen. The walkers are propa- 
gated from X to X' according to the kernel Gx^x'{t) 
in the form (|43[) . In the new sample point X', the ex- 
pectation values are calculated and the starting value of 
Et is adjusted. Eventually, the ground state expectation 
values are projected and sampled. 

The problem is that the wave function ^(X) obey- 
ing the many-body equation (|23|) is complex because the 
Hamiltonian (|22|) contains the magnetic flux. The prob- 
lem can be partly eliminated by means of the fixed-phase 
approximation . Using 



*(X) = |^(X)| e i<l5(x) 



(44) 



with Hamiltonian 



fi 2 N ri 1 

1=1 1 



(39) 



where X = (xi,X2, ■ ■ ■ ,%n), the potential energy V(X) 
incorporates all single-electron interactions and all pair 
electron-electron interactions, and the wave-function 
^(X) is assumed to obey the cyclic condition (|24|) . If 
^(X) is a real function, the DMC method is capable to 
find the exact ground-state solution of equation (|38p . We 
briefly summarize how the DMC worksi^. 

Instead of directly solving the equation (|38j). the DMC 
solves the time-dependent diffusion problem^ 6 - 



(40) 



dV(X,t)/dt= [H-E T )$(X,t), 



where Et is a trial energy. One can write (|4Tj)) in the 
integral formi& 



tf(X,t + r) = J dX' Gx^x'(r) *(X',t), 



(41) 



where G x ^x<(t) = (X| exp(-(7i - £ t )t)|X') is the 
Green's function and r is a time step. Equation (|4ip 
allows to project the ground state ^o(X) from any trial 
wave function \ty(X) which has a nonzero overlap with 
^o(X). Inserting any trial Wy(X) and Et into l|4Tj) gives 

*t(X, r -> oo ) = *o(X)(* |*r) Km exp[-T(E -E T )]. 

T—KX> 

(42) 

By adjusting Et to equal Eq one can make the exponen- 
tial factor constant. It is also important that the Green's 
function Gx^x'(t) can be expressed 16 for r — > as 



G 



(r) a (2irT)- m / 2 



exp 



{X-X>) 



IY2 



2t 



xexp(-T[V(X) + V(X')-2Er]/2). (43) 



one can split the equation (|23p into the real part and 
imaginary part. The real part reads 



fl«ff|tt(X)| = £|*(X)| 



(45) 



where 
H c ff 



N , n ,x 2 

e<p 



2m ^ dxf 2m ^ \ dx t 

1=1 1 i—1 x 



L 



N N 

+7 J2S(xi) + -J2 V(xi- Xj ), (46) 



i,3=l 



and the imaginary part is 

N 



T- 

^ dx { 



I* 



, ; . d _ ecj) 

h dx~^ + T 

UJUi J-J 



= 0. 



(47) 



Equation (|45|) is the effective Schodinger equation for the 
modul |^(X)|, with the Hamiltonian H c g depending on 
the phase $(X). Since |^(X)| is real and H e s has the 
same form as (|3"9"|) . the equation (|4"ij)) can be solved by 
means of the DMC if the phase $(X) is given. In this 
work we choose the trial wave function ^t(X) to be equal 
to the Slater determinant (f2"9"|) of the self-consistently de- 
termined Hartree-Fock ground state and we fix the phase 
$(X) to the phase of this Slater determinant. The DMC 
thus gives the lowest possible ground-state energy E 
within the chosen phase 17 . Eventually, we obtain the 
persistent current from (|28p by finite differences evalu- 
ated by correlated sampling^. Our preliminary DMC 
results are briefly discussed 

in 26 i 27_ 

To go beyond the fixed phase approximation one has 
to couple the DMC solution of equation (|4"5)l with the 
numerical solution of equation (|47|) . We do not attempt 
to do so. However, in section III we check the fixed phase 
approximation by comparing the DMC results with the 
CI calculations which are free of such approximation. 
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E. Configuration-interaction (CI) model 

Before starting with the CI procedure we need to solve 
the single-electron problem 



2m 



ld_ 

i dx 



2tt 4> 
L fa 



j5(x) 



a , 

o 6 — 

a 5 — 

U 4 

u g 3 — 

I 



1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 

quantum number n of Slater determinant x„ 



(48) 

with the cyclic condition ip n (x + L) = ip n (x). We do so 
by using the method described in Section II. A. We obtain 
the wave functions ip n {x) and energy levels e n not only 
for the N lowest energy levels but also for the infinite 
ladder of excited states, of course, with an upper cutoff. 
Consider now the non-interacting many-body problem 




FIG. 1: Schematic sketch of how the determinants \n are 
created for n = 0, 1, . . . in the three-electron system with the 
ladder of single-particle levels restricted to six lowest levels. 
The occupied levels constituting the state \ n are labelled by 
circles. In the CI method the many-body state ^ of the in- 
teracting system is expanded as * = coXo + C1X1 + C2X2 + • ■ ■ • 



where n = 0,1,..., M. This system determines the 



+ r y6(xj)\ Xn(X-) — £„x«(X), eigenvalues E\ and eigenvectors (cq,c' 



1 



) for I = 



(49) 

where X = (x\, x%, . . . , xjy). Clearly, the eigenenergies 
£ n are given as 

and the wave functions Xn are the Slater determinants 
T/V(xi) ... i> nN (xi) 



1 



Xr. 



tp ni (x N ) ... 1pn N (x N ) 



(51) 



where the quantum numbers n±, . . . , njv label the state of 
the first, . . . , A^-th electron, respectively, and the quan- 
tum number n labels the many-body state correspond- 
ing to the specific set of N occupied single-electron levels 
e ni , . . . , £ nN ■ For instance, the figure [1] depicts creation 
of all Slater determinants in the three-electron system, 
when only six lowest single-electron levels are considered. 

To solve the interacting many-body problem the 
CI metho d 28 ' 29 relies on the expansion 



* = cqXo + C1X1 + C2X2 + 



(52) 



Using this expansion and equation (xnlH^) = (xnlE^) 
we obtain from (|23|) the infinite set of equations 



where 



Ec r , 



N 



0,1, 



,00, 



Vnj = \ {Xn\ V ( Xk ~ X l)\*j)- 



(53) 



(54) 



k,l — 1 



We reduce the infinite number of single-energy levels 
Sj to the finite one by introducing a proper upper energy 



cutoff. This reduces the infinite number of equations 
to a certain finite number M+l. We get the finite system 



M 
3=0 



~\~ Vnj ) Cj — ^ 1 



(55) 



0,1,..., M. We obtain the groundstate energy E[ = q and 
groundstate wave function 4 , ; = o by solving the system 
with ARPACK, or alternatively with LAPACK. Once E 
and ^0 ar e known, the persistent current can be obtained 
from (|28f or (f25|) . 

The problem of the CI method is, that for N max single- 
particle levels the expansion (|52|) still contains ^ Ar ™ x 

Slater determinants. For a reasonably chosen cutoff, the 
solution of the equations (|5o| is not feasible already for 
systems with about ten electrons, as the memory require- 
ments are too large. The question is which determinants 
have to be kept in the expansion (|52[) and which can be 
omitted without damaging the final results. This prob- 
lem is known in the quantum chemistr y 28 ' 29 where the 
CI is applied to calculate the energy spectra of molecules. 
We cannot apply here directly the tricks from the quan- 
tum chemistry because our problem is rather special; we 
search for the persistent current in the ID ring pierced by 
magnetic flux. We prefer the following two approaches. 

Our first CI approach, below referred as the FCI, con- 
ceptually resembles the full CI and can be understood 
by means of the sketch in the figure [T] After choosing 
the cutoff for the single-particle energy we add into the 
expansion (|52[) first the Slater determinant of the ground 
state, then all determinants with a single excited electron 
(the single-excitations), all determinants with two ex- 
cited electrons (the double-excitations), all determinants 
with three excited electrons (the triple-excitations), etc. 
If the results do not change, we cut the expansion by 
stopping to increase the number of the excited electrons. 
Some preliminary FCI results are briefly discussed in 27 . 

Our second CI approach, introduced by two of us in 
Ref^S, is called the bucket-brigade CI method (BBCI). 
After choosing the cutoff for the single-particle energy, 
we need to asses importance of the remaining Slater de- 
terminants in the expansion (|52p. In principle, one can 
compute for each determinant Xn the energy (xn \ H |x n ), 
where H is the Hamiltonian (|22p . The obtained numer- 
ical values (Xn\H\xn) can be ordered increasingly and 
we can assume that this orders the determinants Xn ac- 
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cording to their decreasing importance. We can order the 
determinants in the expansion (|52p according to this im- 
portance criterion (a-priori measure of importance) and 
we can truncate the expansion after reaching convergence 
of the final results. This is the basic idea of the BBCI 
method, the term " bucket brigade" is related to the tech- 
nical implementation^ reviewed in the Appendix C. 

In the next section, reliability of our importance cri- 
terion is supported by a successful convergence of the 
BBCI results. We cannot prove the criterion exactly, but 
we can give a simple intuitive motivation. We search for 
the best minimum of the ground-state energy. Therefore, 
the larger the energy {Xn\H\x n ) the smaller should be 
the weight of the given \ n in the ground-state expansion 
([52")) . It is important that the BBCI allows us to treat 
more electrons than the FCI, as it reduces the expansion 
([52]) more efficiently. Another CI algorithm, more effi- 
cient than the FCI, is the CI approach of Ref. [3l|, where 
the determinants are selected by means of a proper Monte 
Carlo algorithm. 



III. RESULTS 
A. Luttinger-liquid scaling in CI and DMC models 

Our calculations are performed for parameters typical 
of a GaAs ring. We use the electron effective mass m = 
0.067mo and electron density n e = N/L — 5 x 10 7 m -1 . 
The strength Vo and range d of our e-e interaction (0 are 
properly varied to demonstrate various e-e interaction 
effects. In practice, screening by metallic gates can be 
used to vary d while Vq can be varied by varying the 
(finite) cross-section of the ID ring. We study the rings 
containing a strong scatterer with \tk F \ 2 <§; 1, because 
in such case the persistent current reaches asymptotic 
dependence on L for relatively small L. Our conclusions 
hold for any ID system, no matter what material is used. 

In figure [5] the persistent current / in the GaAs ring 
with one strong scatterer is calculated as a function of 
the ring length L for magnetic flux <fi = O.250o. The 
scatterer is the <5 barrier with transmission \tk F \ 2 = 0.03. 
The power law LI cx L~ a decays in the log scale linearly 
with slope —a, the same decay show for large L the CI 
and DMC data. Agreement of the FCI and BBCI data 
is very good, the DMC data are close to the CI data. 

Let us check whether our numerical values of a agree 
with the formula ©■ The Fourier transform of our e-e 
interaction (0) is V(q) = 2V d/(l + q 2 d 2 ). Using this 
expression we can write the RG formula (|U) in the form 



N 



16 



32 



64 



OtRG 



4V mk F d 3 
irh 2 (l + 4k F 2 d 2 )' 



(56) 



The dashed lines in the figure [2] show the power law 
LI cx L~ a , with a given by the formula ^ and cxrg cal- 
culated from the formula (|56p . The proportionality factor 
(specified later on) causes in the log scale only a vertical 
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FIG. 2: Persistent current LI jevF versus the ring length L for 
the GaAs ring with a single scatterer. The transmission of the 
scatterer is |tfe F | 2 = 0.03, magnetic flux is <f> = 0.25</>o- The 
upper horizontal axis shows the electron number N = n e L, 
where n e — 5 x 10 7 m" 1 . The range of the e-e interaction 
is fixed to d = 3nm while the magnitude Vo is varied. The 
full curve is our numerical result for Vo = (non-interacting 
electrons). This curve saturates for large L exactly at the 
value 0.5|tfc F |, predicted by the asymptotic formula |TJ. The 
symbols show our FCI, BBCI and DMC results for various 
Vo shown in table [I] the dotted lines connecting the BBCI 
data are a guide for eye. The dashed lines show the Luttinger 



liquid asymptotics LI oc L a for a = (1 + 2o.rg) 



1/2 



1, 



where the values of ana (listed in the figure) are calculated 
from the formula (|56|l . The input parameters Vo and d and 
the resulting ana and a are summarized in the table [I] 



shift of the linear law log(LJ) = — a \og(L) + const. We 
thus can conclude that the CI and DMC data in the fig- 
ure [2] exhibit the L~ a decay with the power a in good 
agreement with the theoretical formulae © and (f5B|) . 

In figure [3] we compare the persistent currents for two 
very different transmissions |£fc F | 2 in order to demon- 
strate that the power a is universal - independent on the 
choice of tk F . Indeed, the CI and DMC data in figure [3] 



Vo (meV) 


d (nm) 


OtRG 




11 


3 


0.0277 


0.0273 


34 


3 


0.0855 


0.0821 


68 


3 


0.171 


0.158 


102 


3 


0.2565 


0.230 



TABLE I: The input parameters Vo and d used in the CI and 
DMC calculations of the figure[2] Also shown are the resulting 
theoretical values of ana and a. 
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FIG. 3: Persistent current LI /evF versus the ring length L for 
the ring with a single scatterer in the regime LI oc L~ a . Mag- 
netic flux is cf> = O.250o, the electron density is N/L — 5 x 10 7 
m _1 . We compare the results for a scatterer with two various 
transmissions |tfc F | 2 . We also make comparison for various 
e-e interactions chosen so that the RG formula (|56p predicts 
for each considered interaction the same ana- Specifically, 
each set (Vo,d) listed in the figure is chosen so that the RG 
formula predicts chrgCVo, d) = 0.2565. The dashed lines show 



1 and 



the formula LI oc L~ a , where a = (1 + 2ccrg 
the proportionality factor is specified in the text. The CI and 
DMC data are shown by symbols. 



exhibit the same a for both transmissions. 

One has to note that the formulae ([4]) and ([6]) are ro- 
bust against the choice of the e-e interaction in the sense, 
that they give the same a for all e-e interactions with the 
same value of [V(0) — V(2Jcf)]/vf- Obviously, this can 
be the case for many various choices of V (x — x'). 

To see whether such robustness exists in our many- 
body model, in figure[3]we also make comparison for vari- 
ous e-e interactions chosen so that the formula (JHSJ) gives 
for various sets (Vo,d) the value <XRG(Vo,d) = 0.2565. 

Our CI and DMC data indeed show for all considered 

l li 

(Vb, d) the decay LI oc L~ a , where a = (1 + 2anci) —1 
and ana — 0.2565. In summary, our methods seem to 
give the same a for various sets (Vb , d) fulfilling the equa- 
tion o>iig{Vo, d) =const. Thus, our methods seem to con- 
firm the above discussed robustness of the formulae ((3]) 
and ^ against the choice of the e-e interaction. 

However, we will show later on that this robustness 
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FIG. 4: Data from figure [3] plotted as the ratio ((59}, where 
I(tk F ,i) and I(tk F ,2) are the persistent currents for lifc^il 2 = 
0.03 and |tfc F ,2| 2 = 0.003, respectively. 



in fact holds only for the e-e interactions which are very 
short-ranged [all sets (Vo,d) in figure [3] actually belong 
to this special limit]. We know that we are in this limit 
as long as our methods confirm the formulae ^ and © . 
This is still the case in the following discussion. 

We specify the proportionality factor in the formula 
LI oc L~ a . Following the Introduction, we replace in the 
single-particle formula (P) the bare amplitude tk F by the 
renormalized amplitude tk F oc (tk F /\fk F \)L~ a . We get 



n ev F \t kF \ r _ Q . , 
k ~^-r 7^ — r L sin (2n 
2L r kF y 



a > 0, 



(57) 



where C is the proportionality factor. Since a is universal 
(independent on tk F ) , the formula (|57[) implies that 



I(tk F ,l) _ \tk f 



\ r k r 



I(t kp ,2) \tkr,2Wkp,2\ 



(58) 



as long as C does not depend on tk F - To see that the 
CI and DMC data in figure [3] fulfill the equality f55]). we 
show these data again in figure [4] in terms of the ratio 



|r fef ,, 2 | I{h F ,i) 



|rfc F ,i| I{tk F ,2Y 



(59) 



which should equal unity for L — > oo. One sees that it in- 
deed approaches unity for the CI as well as DMC results. 
More precisely, the CI data show a few percent deviation 
from unity which tends to disappear with increasing L. 
A small deviation from unity show also the DMC data 
where convergence towards unity is not clear due to the 
stochastic noise of the Monte Carlo method. 

Thus, if we ignore a small finite-size effect in figure 21 
our data in figure [3] are in accord with equation 
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FIG. 5: Dependence Co (a) obtained by fitting the formula 
(1601) to our CI and DMC data for Ll/eVf. Symbols show the 
fitted values of Co- The dashed line is the function Co (a) = 
1 + 1.66a which fits the presented data points. 



Therefore, we can fit our data in the L~ a regime by the 
asymptotic formula (|57|) . For the purpose of fitting it is 
useful to define C =n e ~ a Co. Setting this definition into 
the formula ([57)1 we obtain 

— = lc (a) M sin(27T0/0o), (60) 
ev F 2 \r kF | 

where Co depends solely on the parameter a (see below) . 
Instead of ([57)1 we use (ftJOl and we fit Co instead of C. 

The dashed lines in figures [2] and [3] show the depen- 
dence fitted to the BBCI data in these figures. 
Specifically, Co is fitted and the function N~ a in ([50)) 

1/2 

is evaluated by using a = (1 + 2a rc) — 1, where 
£*RG(Vo,d) is given by the formula ([56"]) . The resulting 
values of Co are shown by open symbols in the figure [5] 
The full symbols in the figure [5] show the values of Co 
obtained when we fit by means of (|60"f the DMC data 
in figures [2] and [31 Also shown are the results obtained 
in the same way for another values of a. Note that for 
all (Vo, d) giving the same a(Vo, d) the obtained values of 
Co are essentially the same. This means that Co depends 
solely on the parameter a. We also note that Co does not 
depend on tk F , as has been documented in the figure [4] 
The linear fit in the figure [SJ 



C (a) = 1 + 1.66a, 



(61) 



should be viewed as a first estimate of the analytical de- 
pendence. To extract a very precise formula for Co (a), 
it would be desirable to simulate larger systems in order 
to better suppress the finite size effect in figure |U 

We conclude that the right-hand side of (|rjD)l depends 
on a single parameter a via the function Co(a)N~ a , 
which is single- valued for various (Vo, d) giving the same 
a(Vo, d). Our results confirm the decay N~~ a , where a is 
given by the formulae © and (|3]). 

However, the formulae ([6]) and ([4]) in fact hold only for 
very small d. In the figure [5] we compare the BBCI data 
from the figure [3J with the BBCI data obtained for an- 
other two sets (Vo, d) with d as large as 12nm and 24nm. 
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FIG. 6: The circles, squares and diamonds are the BBCI data 
from the figure [3] the triangles and inverted triangles are the 
BBCI data for another two sets (Vo,d). For all considered 
(Vo,d), the formula (I51)[) gives aRo(Vo,d) — 0.2565 and the 
formula © predicts a = (1 + 2aRc) 1 ^' 2 — 1 = 0.23. However, 
the dependence LI oc L -0 ' 23 (the dashed curves with different 
offset) can be seen to fit only the BBCI data for d < 3nm. The 
BBCI data for d — 12nm are excellently fitted by the curve 
LI oc L -0 ' 222 [the full line] and the BBCI data for d = 24nm 
are fitted by LI oc £~ ' 207 [the dotted line]. This shows that 
the formulae (|4]) and ([6]) hold only for the e-e interactions 
which are very short-ranged (see the text). 



All considered (Vo,d) are chosen so that the formulae 
([5*6]) and (0 give for each (Vo,d) the same value of a, 
in particular a = 0.23. However, it can be seen, that 
the dependence LI cx L~ 23 fits only the BBCI data for 
d < 3nm. The BBCI data for d = 12nm and d = 24nm 
are excellently fitted by the curves LI cx L~ - 222 and 
LI cx L~ - 207 , respectively, i.e., with increasing d our 
numerically obtained a decreases albeit the value of 
ajic(Vo, d) is fixed. This is a strong indication that the 
formulae ([5]) and ^ hold only for small d. We support 
this numerical finding by the following analytical proof. 

In the appendix D we prove analytically that the ma- 
trix elements of the e-e interaction are independent on d 
for d — * at a fixed value of chro^o, d). In other words, 
for d — > the matrix elements depend on a single pa- 
rameter olrq, otherwise they depend on two parameters 
Vo and d. According to the appendix, the limit d — > 
means 4k pd 2 <C 1. In our calculations 1/2 k p — 3nm and 
the currents in the figure [B] are close to each other for all 
d < 3nm. If we increase ana (for instance as in the table 
Inl and figure [JJ, we need to take d < l/2kp to see the 
rf-independent current convincingly. 

Note also another result in figure [BJ As d exceeds 
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Vo (meV) 


d (nm) 


CtRG 


a 


symbol in Fig. [7] 


11 


3 


0.0277 


0.0273 


o 


34 


3 


0.0855 


0.0821 


o 


68 


3 


0.171 


0.158 


o 


1068 


1 


0.171 


0.158 


□ 


102 


3 


0.2565 


0.230 


o 


1602 


1 


0.2565 


0.230 


□ 


11957 


0.5 


0.2565 


0.230 


A 


15942 


1 


0.342 


0.298 


□ 


15942 


0.5 


0.342 


0.298 


A 


3142 


1 


0.5 


0.414 


□ 


23445 


0.5 


0.5 


0.414 


A 



TABLE II: The input parameters Vo and d used in the BBCI 
calculations of figure [7] and the values of ckrg and a resulting 
from the formulae (|56[1 and (JB]). The last column of the table 
ascribes a symbol to each set (Vo,d, ocrg)- These symbols are 
used in figure [7] to show the BBCI data. 



l/2kp the persistent current still decays like LI «L _ ", 
but a depends on two parameters Vq and d rather than 
on a single one, ana- This regime is not captured by the 
formulae ((6]) and ((4|. We have obtained similar results 
(not shown) also for other values of ana ■ 

Consider now a broader range of the e-e interaction 
parameters Vo and d, listed in the table ITT1 The resulting 
persistent currents for these parameters are shown in the 
figure [7] All shown BBCI data asymptotically converge 
to the Luttinger-liquid behavior described by the formula 
([6"0")) . with the power a in accord with the formulae © 
and Q. We wish to stress the following aspects. 



We can write (l60l) as 



2LI 



jX. 



where x = N a 



ev F C \r k p\' 

and Co (a) is given by (|6T|) . Similarly, the BBCI data in 
the asymptotic regime can be normalized as 2LI /evpCo 
and plotted as a function of the variable x = N~ a . This 
is done in inset to the figure [7J Indeed, all BBCI data 



That 



in inset collapse to a single curve f{x) 

a single linear curve involves the BBCI data for many 
various N, Vo and d, is a clear sign of the Luttinger liquid 
with e-e interaction depending on one parameter q.rg- 
It also documents that our calculations are reliable for a 
broad range of variables N, Vo and d, the function Co (a) 
is certainly determined with reasonable accuracy. 

An interesting finding in figure [Jj is that for large a^a 
the BBCI data show the LI oc L~ a decay already for ten 
electrons. In other words, the Luttinger-liquid behavior 
arises in the system with only ten particles. For compar- 
ison, in the lattice model o 12 ' 13 i 14 ' 15 the asymptotic power 
law is observed for a much larger number of electrons. 

If we compare the formulae (frjO"]) and (fT]), we see that 
the interaction modifies the amplitude tk F as 



t kF c {t kF /\f kF \)N' a = (h F /\f kp \){n- x /L) a . (62) 
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FIG. 7: Persistent current versus the ring length for the same 
ring as in the figure [2] but for a broader range of the interac- 
tion parameters Vo and d (table [II]). The circles, squares and 
triangles show the results of the BBCI method for various sets 
(Vo, d, org) listed in the table [TT1 The dotted lines are a guide 
for eye. The dashed lines show the asymptotic law LI oc L~ a 
plotted in the form (|60|) . where a = (1 + 2o:rg) 1 ^ 2 — 1 and 
cxrg is obtained from ([56} ■ The BBCI data in the asymptotic 
regime are selected, normalized as 2LI / evFCo [with Co(a) 
given by the formula (|61[) ] , and plotted in inset as a function 
of the variable x = N~ a . The full line is the dependence 

f(x) = |-^ F | x, predicted by the formula (|60[) . 



This result scales with the length nj 1 while the RG re- 
sult ([3|) scales with d. What is the origin of this differ- 
ence? The RG result ([3]) should hold for any d obeying 
the inequality ([S]), but the inequality (O can in princi- 
ple be fulfilled also for d < l/2fep, i.e., just for those d 
for which we have obtained the formula ([62]) . We think 
that the difference between the two results is due to the 
different physical models. The formula (p)2")) holds in our 
continuous model, with the single-particle energy disper- 
sion being truly parabolic. The RG formula holds 
in the model^, with the energy dispersion linearized in 
the interval < Sf — hvp/d, Ef + hvp /d > around the 
Fermi energy. For small d the band width hvp/d be- 
comes comparable with ep and the model^ has a truly 
linear dispersion, like the Luttinger-liquid model^ which 
also shows scaling by factor (d/L) a . For d » l/2kp 
the linearization is unessential because of hvp/d <C 6p, 
hence both models should give the same results. The 
limit d l/2kp is however not feasible by our numeri- 
cal methods for computational reasons. 
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FIG. 9: The Hartree-Fock data from the figure [8] plotted as 
the ratio ()59|) . where I(tk F ,i) and I(tk F ,i) are the persistent 
currents for the transmission probabilities |i fcj^ , i | 2 = 0.03 and 



\tk f 



0.003, respectively. 



FIG. 8: Persistent current LI /evF versus ring length L for the 
ring with single scatterer, calculated in the Hartree-Fock ap- 
proximation for parameters (f> = 0.25</>o, N/L = 5 x 10 7 m -1 , 
d — 3nm, and various Vb and ifc F | 2 as listed in the figure. 
The left panel shows the nonlocal Hartree-Fock results, the 
right panel shows the Hartree-Fock results in the local Fock 
approximation (see the text). The Hartree-Fock results are 
shown by symbols, the dashed lines show the Luttinger liquid 
asymptotics LI oc L~ a , where a = (1 + 2q_rg) 1,/2 — 1 and 
a RO is given by the RG formula ||SB}: a RG = 0.0277, 0.0855, 
and 0.1710 for Vo = 11,34, and 68 meV, respectively. 

B. Universal scaling in the Hartree-Fock model 

Now we discuss the persistent currents obtained in the 
self-consistent Hartree-Fock approximation (Sect.IIC) 
which ignores correlations except for the Fock exchange. 

In figure [8] we show the I(L) dependence, calculated 
in the self-consistent Hartree-Fock approximation for 
the ring parameters already encountered in our corre- 
lated many-body calculations. The left panel shows 
the Hartree-Fock data obtained for the Fock term ([33]) 
treated as is, the right panel shows the Hartree-Fock data 
for the Fock term (|33|) approximated a a 32 ' 33 

U F (x)~-^2 J dx'V(x - x')Re [ij}*,(x')il> n >(x)\ ■ (63) 

n' 

Unlike the nonlocal interaction ()33p the approximation 
(|6"3")) is local - it does not depend on n. One readily ob- 
tains ([63)) from (|33|) by applying the 'almost closure re- 
lation' Yin' ^n'i^^n'ix) — 5{x — x'). Once this approx- 
imation is adopted, the final expression for the current 
has to be approximated correspondingly 34 . 

Let us compare our Hartree-Fock calculations with the 
Luttinger liquid asymptotics LI oc L~ a . As can be seen 
in the figure El both Hartree-Fock approaches show the 
I(L) dependence decaying for large L faster than L~ 1-a . 



We will see below that for L — > oo the decay is in fact ex- 
ponential. Only for the weak e-e interaction the Hartree- 
Fock data show the decay L _1_a , reported also in our 
previous Hartree-Fock studie o 35 i 36 . However, also in this 
case we expect exponential decay at very large L. 

Now we show that also in the Hartree-Fock approxima- 
tion the slope of log(7(L)) becomes for L — > oo universal 
- dependent only on the e-e interaction. Following equa- 
tions ([571) - (|59l . such universality means that the ratio 
\tk F ,2\/\rk F ,2\ i(tk F ,i) a pp roac h es um ty for L — > oo. If we 

\tk F ,l\/\r kF ,i\ I{tk F , 2 ) J 

redraw in terms of this ratio the Hartree-Fock data from 
figure El we see (in the figure [5]) , that the ratio converges 
with increasing L to unity. The stronger the interaction 
the better the convergence in the figure [9l for further im- 
provement a further increase of L is needed. Universality 
of the Hartree-Fock results in figure [5] is quite similar to 
the universality of the correlated results in figure |U 

Moreover, we would like to show that the Hartree-Fock 
model resembles the correlated model also in the follow- 
ing respect. In the limit d — > the I(L) curve is de- 
termined by a single e-e interaction parameter chrc [i.e., 
I{L) is the same for various (Vq, d) giving the same value 
of anc(Vo, d)]. This is demonstrated in the figure ITTJl 

Figure [TO] shows the Hartree-Fock results for the weak 
and strong e-e interactions. Also shown are the corre- 
sponding correlated results from the figure [H in partic- 
ular the BBCI data (full lines) and the Luttinger-liquid 
curves (dashed lines). Note the following details. 

For weak e-e interaction (ana = 0.0277) the Hartree- 
Fock data are robust against various (Vo,d) obeying the 
equation cenG(Vo,d) = 0.0277. This is illustrated by the 
agreement of the Hartree-Fock data for (Vb = 173meV, 
d = lnm) and (Vb = llmeV, d = 3nm). Note also that 
these Hartree-Fock data almost reproduce the correlated 
results at sizes L considered in the figure. To obtain 
conclusions valid for any L and any interaction strength, 



12 



10" 



S 10' 



N 

10 100 | 

non-local HF 

gS°<^) = °-°277 




AV = 1 1956 meV, rf = 0.5 nm 
V = 173 meV, d = 1 nm 
V = 1603 meV, d= 1 nm 

• V = 11 meV, d = 3 nm 



o - 



0.02 



102 meV, d = 3 nm 



0.2 2 
L (um) 



AT 

10 100 

local HF 

fe^W) = o.02 77 




1 A 



>2. 

a 



0.02 



0.2 2 
L (urn) 



FIG. 10: Length-dependence of the persistent current in the 
ring with single scatterer, obtained in the Hartree-Fock ap- 
proximation. The Hartree-Fock results are shown by sym- 
bols. The ring parameters are: <j> = 0.25<?io, N/L = 5 x 10 7 
m _1 , and |~fc F | 2 = 0.03. The weak and strong e-e interac- 
tions are simulated for various (Vo,d) obeying the equations 
oma{Va,d) = 0.0277 and aRc(Vo,d) — 0.2565, respectively, 
where a. no (Vq, d) is given by the RG formula (|56|l . Also shown 
are the corresponding BBCI data (full lines) and asymptotic 
LI oc L~ a curves (dashed lines), taken from the figure [2] 



we have to look at the stronger e-e interaction. 

For strong e-e interaction (a = 0.2565) the Hartree- 
Fock data fail to agree with the correlated results. How- 
ever, the Hartree-Fock data still tend to be robust against 
various (Vq, d) obeying the equation a.RG(Vb, d) — 0.2565 
in the limit d — > 0. This universal dependence on a single 
e-e interaction parameter ana is clearly visible both in 
the nonlocal and local Hartree-Fock calculation. 

Finally, we would like to show that the decay of I(L) in 
the Hartree-Fock approximation is exponential for large 
L. This exponential decay is demonstrated in the figure 
lilt where the I(L) dependence is presented in the semi- 
logarithmic scale. We need to consider large L and strong 
e-e interaction, which strongly prolongs the time of the 
Hartree-Fock computation. Fortunately, the calculation 
is feasible in the local approximation ([55]) . 

The I(L) curves in figure fTTI do not decay equally fast, 
albeit the value of ctnc(Vo,d) is the same for all consid- 
ered (Vq, d). It can however be seen that if we keep the 
same value of aRG(Vo,d) in the limit d —> 0, then the 
slope of I{L) is determined solely by the value of oirg- 
This is again a clear manifestation of the universal de- 
pendence on a single e-e interaction parameter. The limit 

— ^ is well represented by the data for d = O.lnm. 

In summary, our self-consistent Hartree-Fock results 
show, that the slope of \og(I(L)) in the limit of large L is 
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FIG. 11: Length-dependence of the persistent current in the 
ring with single scatterer, obtained in the local Hartree-Fock 
approximation. The Hartree-Fock data are shown by symbols. 
The ring parameters are: <j) = O.250O, N/L = 5 x 10 7 m -1 , 
and |~fc F | 2 = 0.03. The parameters of the e-e interaction, Vb 
and d, are listed in the figure. The considered sets (Vo,d) 
obey the equation aRc(Vo,d) = 0.2565, where a.Ra{Vo,d) is 
given by the RG formula (|56[) . For large L the current decays 
with L linearly, but note that the scale is semi-logarithmic. 
The dotted lines fit linearly the Hartree-Fock data at large L. 



still universal - dependent solely on the e-e interaction not 
on the strength of the scatterer. Moreover, for d — > this 
dependence on the e-e interaction is determined by a sin- 
gle e-e interaction parameter, specifically by otRG given 
by the RG formula (|56|) . These features are very similar 
to the universal behavior observed in our correlated cal- 
culations. A major difference is that in the Hartree-Fock 
approximation the asymptotic decay of I(L) is exponen- 
tial with L, as we have shown numerically. 



C. Reliability and limitations of the CI and DMC 

Both DMC and CI methods in their most rigorous for- 
mulations show exponential scaling of the computer time 
in the number of particles. In CI this is due the ne- 
cessity of infinite size of the basis and high level of ex- 
citations which are required to make the method size 
consisten t 28 ' 29 . In practice, for not too large systems, 
useful and accurate results can be obtained providing a 
sufficiently large set of virtual orbitals can be treated. 
The question is whether the low-order (doubles, triples, 
quadruples) excitations are enough to capture the key 
many-body effects. For the DMC method, the exponen- 
tial scaling originates in the fermion sign problem which 
is in practice avoided by the fixed-node or fixed-phase 
approximations. The usefulness of the method crucially 
relies on accuracy of the trial functions and ability to 
efficiently describe the impact of the many-body effects 
on accuracy of the nodes or phase 16 . The DMC practice 



0.11 




FIG. 12: Persistent current versus magnetic flux for the GaAs 
ring with N = 4 and L — N/n e = 80 nm. The transmission 
of the scatterer is \tk F | 2 = 0.03. The e-e interaction range is 
set to d = 3 nm and the e-e interaction magnitude is Vo = 
68 meV in the left panel and Vo — 200 meV in the right 
panel. The triangles show the results of the fixed-phase DMC 
calculation with the trial wave function ^t(X) equal to the 
Slater determinant ([29} of the self-consistently determined 
Hartree-Fock ground state. The squares show the results of 
the fixed-phase DMC calculation with the trial wave function 
$t(X) equal to the Slater determinant of the non-interacting 
ground state. The CI results are shown by circles connected 
by the dotted line. For completeness, the results for the non- 
interacting ground state are shown in a full line and the results 
of the self-consistent Hartree-Fock calculation are shown in a 
dashed line. 

shows that very large systems can be treated neverthe- 
less the fixed- node/phase bias is present and it depends 
on phenomena of interest whether it can affect the re- 
sults. 

In subsection IIIA, the largest many-body systems (48 
electrons) were simulated by the CI method while only 
at most 32 electrons were simulated by the DMC. This 
is caused by the fact that the phase of the Hartree- 
Fock wave function becomes less accurate with increasing 
strength of the interaction and also with increasing size of 
the system. This is not too difficult to understand since 
in larger and strongly interacting systems the collective 
excitations become more complicated and the single de- 
terminant wave function becomes very poor representa- 
tion of the actual ground state. The many-body effects 
description built into the trial function thus directly de- 
termines the accuracy of the obtained results. 

In figure[T2]we show the persistent current versus mag- 
netic flux for the ring with four electrons. The range of 
the e-e interaction is d — 3nm, the e-e interaction magni- 
tude is Vq = 68mcV in the left panel and Vo = 200meV 
in the right panel. The triangles represent the DMC cal- 
culation with the trial wave function \I/t(X.) equal to 
the Slater determinant (f2U]) of the Hartree-Fock ground 
state. The squares show the DMC calculation with the 
trial wave function equal to the Slater determinant of 
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FIG. 13: Comparison of various phases from the DMC 
calculations of figure [12] at (j> = 0.2500- The top panel 
shows the phase $(xi, X2, xz, Xa) of the non-interacting trial 
wave function ^t(xi, X2, xz, xa) evaluated at xz — —L/6 
and xa = L/6. It is labelled as $(0meV) since it holds 
for Vo — OmeV. The middle panel shows the phase differ- 
ence |$(68meV) — $(0meV)|, where <E>(68meV) is the phase 
$(a;i, X2, xz = —L/6,X4 = L/6) of the trial wave function 
< &t(xi,X2,xz,xa) obtained for Vo = 68meV by the Hartree- 
Fock method. The bottom panel shows the phase difference 
|$(200meV) — $(0meV)|. We note that the abrupt change of 
the phase for Xi — Xj arises because the Slater determinant 
changes abruptly its sign for Xi «-» Xj. 



the non-interacting ground state. The results of both 
DMC calculations are very close when Vo — 68meV, but 
they are very different when Vo = 200meV. This suggests 
that the DMC results for Vo = 200meV are strongly af- 
fected by the phase of the trial wave function. Indeed, 
for Vo = 68meV both DMC calculations essentially agree 
with the corresponding CI calculation (circles), but for 
Vo = 200meV the differences with the CI data are very 
clear. 

It is interesting to see the origin of the fixed-phase 
biases. For illustration, in figure [T2] we compare the 
phases involved in the DMC calculations of figure [T2J 
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FIG. 14: Convergence of the persistent current in the FCI 
calculation in dependence on the number of the considered 
single-electron levels for the single-excitations, single- and 
double-excitations, etc. The top panel shows the currents 
obtained from the formula / = {^ol / l^o), the bottom panel 
shows the currents obtained in the same FCI calculation from 
the formula / = —dEo/d<f>. The results are shown by sym- 
bols. They were obtained for the ring with parameters N = 8 
and L = N/n e = 0.16 fj,m, penetrated by magnetic flux 
4> = 0.2500- The transmission of the scatterer is \tk F | 2 = 0.03. 
The e-e interaction is given by Vo = 102 meV and d = 3 nm. 



The phase difference |$(68meV) — <£>(0meV)| is small and 
this is why the corresponding DMC currents in the left 
panel of figure fT2l almost coincide. The phase difference 
|$(200meV) — <I>(0meV)| is large and the corresponding 
DMC currents in the right panel of figure [T2l clearly differ. 

We conclude that the main limitation of our DMC re- 
sults is the Hartree-Fock approximation of the correct 
phase, which deteriorates with increase of the e-e inter- 
action and system size. In particular, we see in the figure 
[T2l that the DMC with the Hartree-Fock trial wave func- 
tion underestimates the persistent current. For instance, 
in figure [2] the DMC data for the ana = 0.2565 exhibit 
a weak underestimation which is of the same origin. 

Reliability of the CI results depends on how the re- 
sults converge with increasing the number of the Slater 
determinants in the expansion 
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FIG. 15: Convergence of the persistent current in the BBCI 
calculation for various N max in dependence on the number of 
the Slater determinants involved in the expansion (|52[) . We 
recall that the determinants are selected by using the bucket- 
brigade algorithm of Appendix C and N max is the number of 
the considered single-electron levels. In this BBCI calculation 
the current is evaluated by using the formula / = { M/o 1 1 l^o) 
and the parameters of the ring are the same as in the FCI 
calculation of the preceding figure. The value 0.0765 is the 
average from the (well saturated) data for N max > 24. 



The figure [TH shows a typical convergence process in 
the FCI calculation. We recall that we add into the ex- 
pansion (|52p first the Slater determinant of the ground 
state, then all determinants with a single excited elec- 
tron (the single-excitations), all determinants with two 
excited electrons (the double-excitations), etc. The fig- 
ure 1141 demonstrates how the persistent current satu- 
rates with increasing the size of the single-electron ba- 
sis (the number of the considered single-electron lev- 
els) for the single-excitations only, for the single- and 
double-excitations, etc. Moreover, the top and bottom 
panel compare convergence of the persistent currents 
calculated by means of the formulae / = (^ol^l^o) 
and I = —dEo/d(f>, respectively. It can be seen that 
both approaches converge to the same final result when 
the single-, double-, triple-, and quadruple-excitations 
are considered. However, precision of the results in 
the bottom panel is good already for the single- and 
double-excitations while in the top panel also the triple- 
excitations are needed to achieve a comparable precision. 

The FCI calculations with the formula I = —dEo/dcj) 
therefore consume much less computer time and mem- 
ory. The FCI data presented in this text were obtained 
by means of the formula / = —dEo/d<j) and by consider- 
ing the single- and double-excitations. In a few cases, 
sufficiency of the achieved convergence was confirmed 
by adding the triple-excitations or even the quadruple- 
excitations, with a similar success as in the figure [T4l 

The figure [TH] shows typical convergence of the BBCI 
calculation in dependence on the number of the Slater de- 
terminants involved in the expansion (|52p . The persistent 
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FIG. 16: The left panel shows how the ground-state energy Eo 
converges in the FCI calculation of figure [T4l The right panel 
shows how Eo converges in the BBCI calculation of figure [151 



currents in the BBCI converge to the value 0.0765, in ac- 
cord with the value reached by the FCI data in figure [T4l 
The BBCI data in this text were obtained by using the 
formula / = (* | I |*o)- The BBCI relying on the for- 
mula / = —dEo/dcf) gives the same results (not shown), 
but the computational time is about twice longer. 

In figure[I2]we demonstrate convergence of the ground- 
state energy Eo in the FCI (left panel) and BBCI calcula- 
tion (right panel) . In both cases Eq converges to a certain 
minimum. Here Eq is better minimized in the FCI calcu- 
lation than in the BBCI, but we note that the BBCI re- 
sult would eventually converge (for a significantly larger 
number of determinants) to the same minimum. In both 
cases, however, we need to add more determinants to 
achieve a perfectly saturated Eq. In contrast to this, 
if we look at the convergence of the corresponding per- 
sistent currents ( figures [T4l and [T5l respectively), we see 
that it is satisfactory finished both in the FCI and BBCI. 
Thus, as long as we are interested in the persistent cur- 
rent, it is economical to look at the convergence of the 
current rather than at the convergence of Eq. 

In figure Q~7] we show the ground-state energy and per- 
sistent current as functions of magnetic flux for the same 
parameters as in figures Q3J [T5j and [THJ The CI and 
DMC results for the currents are in good agreement. 
However, the CI and DMC ground-state energies exhibit 
small differences which deserve a comment. In particu- 
lar, in spite of the fixed-phase approximation the DMC 
ground-state energy shows the best minimum. The rea- 
sons for this are the following. First, the FCI results in 
the figure [T71 were obtained by including the single, dou- 
ble, and triple-excitations (see the discussion of figures 
[T4l and [16]) . The quadruple-excitations, etc., would shift 
the FCI ground-state energy to a slightly lower value. 
Second, the presented BBCI results were obtained for 
Nmax = 64 and for 1.6 x 10 4 Slater determinants (see fig- 
ures [T5] and . Inclusion of more determinants would 
shift the BBCI ground-state energy to a lower value. 

It is tempting to think that the result showing the low- 
est ground-state energy is the best one. Such criterion 
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FIG. 17: The ground-state energy Eo and persistent current 
as functions of magnetic flux for the ring with parameters N = 
8 and L = N/n e = 0.16 (j,m. The transmission of the scatterer 
is \tk F \ 2 = 0.03. The parameters of the e-e interaction are 
Vo — 102 meV and d = 3 nm, they give ana = 0.2565. 
The CI, DMC and Hartree-Fock (HF) results are shown by 
symbols. The full lines are a guide for eye, the dashed line 
shows the formula (|60[l . with a = (1 + 2anc) 1 ^ 2 — 1 in the 
function N~ a and with the value of Co taken from figure 



indeed follows from the variational principle if we com- 
pare various methods applied to the same Hamiltonian. 
Here the CI methods are applied to the exact Hamilto- 
nian but the DMC is applied to the effective Hamilto- 
nian. If the ground-state energy of the effective Hamil- 
tonian is lower than the ground-state energy of the exact 
Hamiltonian, this is not in conflict with the variational 
principle. It should also be stressed that the conver- 
gence of other properties such as expectation value of 
the current operator is not necessarily governed by the 
same behavior as the energy convergence. We believe 
that our CI results are closer to the true ones due to 
the good CI convergence demonstrated above. On the 
other hand the crude Hartree-Fock wave functions used in 
the fixed-phase approximation are probably not accurate 
enough to provide accurate currents for strongly inter- 
acting limit. This remains an interesting point for future 
studies since more accurate trial wave functions can be 
constructed by means of pairing orbitals and pfaffians, 
expansions in pfaffians and/or using backflow (dressed) 
many-body coordinates^. 



IV. SUMMARY AND CONCLUDING 
REMARKS 

We have studied the persistent current of the inter- 
acting spinless electrons in a continuous ID ring with 
a single strongly reflecting scatterer. We have included 
correlations by solving the Schrodinger equation for sev- 
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eral tens of electrons interacting via the pair interaction 
V{x — x') — Vq exp(— \x — x'\ /d). Our aim was to solve 
this continuous many-body problem without any phys- 
ical approximation and to examine microscopically the 
power-law behavior of the Luttinger liquid. We have 
used advanced CI and DMC methods which, unlike the 
Luttinger-liquid model, do not rely on the Bozonization 
technique. 

In the past, similar studies were performed by numeri- 
cal RG methods in the lattice mode l 12 ' 13 ! 14 ' 15 . Our meth- 
ods do not rely on the RG techniques and thus serve as 
independent confirmation of such approaches. Moreover, 
dealing with the continuous model, we can vary the range 
of the e-e interaction in order to test the robustness of 
the Luttinger-liquid power laws against various shapes of 
V(x — x 1 ). This point was not addressed in the lattice- 
model studies as the range of the interaction was usually 
fixed to the on-site interaction and/or to the nearest- 
neighbor-site interaction. Our major findings are: 

(i) Our CI and DMC calculations confirm that the per- 
sistent current exhibits the asymptotic power-law behav- 
ior of the Luttinger liquid, / oc L -1- ". 

(ii) In our continuous model, the question whether the 
power a is determined by a single e-e interaction parame- 
ter ana is addressed by using various shapes of V(x — x') 
giving the same value of o.rg- Our numerical values of 

1/2 

a confirm the theoretical formula a = (1 + 2a rg) ~ 1 
with ana given by the RG expression (U) , but only if the 
range of V(x — x 1 ) is small (d < l/2fcp). 

(iii) The CI data for a>RG ^ 0.3 show onset of the 
asymptotics / oc L _1 ~ Q already for ten electrons. In 
other words, the Luttinger-liquid behavior emerges in the 
system with only ten particles. For comparison, a far 
much larger number of electrons is needed to observe the 
asymptotic power law in the lattice mode l 12 ' 13 ' 14 ' 15 . To 
understand origin of this difference, it would be desirable 
to study smooth transition from the lattice model to the 
continuous model. 

(iv) We have treated the e-e interaction in the self- 
consistent Hartree-Fock approximation. We observe for 
large L the exponentially decaying I(L) instead of the 
power law. However, the slope of log(/(L)) still depends 
solely on the parameter olrg given by the RG formula 
dH), as long as the range of V(x — x') approaches zero. 

(v) We have discussed the reliability and limits of our 
CI and DMC calculations. The CI methods appear to 
provide convergent results for the sizes and strengths of 
interactions we have studied. The reliability of the FCI 
and BBCI results is easy to analyze in this case and both 
methods converge to the same results although the BBCI 
can treat larger systems. The main limitation of our 
DMC calculations is the accuracy of Hartree-Fock trial 
function phase which has been employed in the fixed- 
phase calculations and appears to be a rather poor ap- 
proximation for strong e-e interaction. 

Although our work is not aimed to address experimen- 
tal aspects, nevertheless, one of our findings might be a 
motivation for experimental work. It might be interest- 
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FIG. 18: Electron-electron interaction versus distance in the 
GaAs ID system. The dotted line shows the bare e-e in- 
teraction Vbarc(£ - x ) = \ x - x i\+ xa for 6 = 12-5e and 
xo = 3 nm, where the cutoff xo mimics the finite wire thick- 
ness. The thin full line shows the potential of a single elec- 
tron screened by the free ID electron gas, calculated in the 
Hartree picture^. The full line is the best fit by formula 
V(x — x') = Voe~\ x ~ x V d , where d = 3nm essentially co- 
incides with the value of l/2kp. We expect that external 
gates would further increase the screening, i.e., it is indeed 
meaningful to focus on the range d < 1/2&F as we did in 
our work. The Friedel oscillations in the figure are artefact 
of static screening and would be further suppressed by the 
gates. 

ing to observe onset of the Luttinger-liquid behavior in 
the ID system with a small number of strongly interact- 
ing electrons, e.g. with only ten electrons as predicts our 
work. In this respect we mention, that the screened e-e 
interaction (|7|) with the parameters Vq and d considered 
in our calculations is quite realistic. This is documented 
in figure IT51 where the model interaction ([7]) is compared 
with the microscopic Hartree screening of the bare e-e 
interaction. 
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Appendix A: Numerical algorithm for calculation of 
amplitudes t k and 

Consider again the segment (-L/2, L/2) embedded be- 
tween two perfect semi-infinite leads. Outside the seg- 
ment the electron wave function reads 

<p k (x) = e lkx + r k e~ lkx , Vk {x) = t k e* kx (64) 
for the electron impinging the segment from the left and 

<P- k (x) = t' h e- ik \ v-k{x) = e~ ikx + r' k e* kx (65) 
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for the electron impinging the segment from the right, 
where k > 0, is the reflection amplitude, and 
is the transmission amplitude. Therefore, the bound- 
ary conditions for elastic tunneling through the segment 
(— L/2, L/2) have a standard form 



-ik4 



rke 



<p k (-)=t k e ik -*, (66) 



¥>-fc(— S") 



¥>-*(ir) 



r / zfc- 



(67) 



To calculate £fc and r^, we have to solve equation pip 
as a tunneling problem with boundary conditions (|66[) . 
We write equation pip in the discrete form (f2"Tj) . Con- 
sider now the electron with fc > 0. It leaves the segment 
(-L/2,L/2) at x = L/2 as a free wave cx e lkL/2 . We 
can thus initialize the scheme (|2ip at x = L/2 + A and 
x = L/2 by using 



, num 



(L/2 + A) = e * fc (i/2+A) ; ^num ^ ^ = e *fcL/2_ 



By means of (|2ip we generate the numerical solution 
(p™ m (xj) starting at x — L/2 — A and ending at £ = 
-L/2. Since y k (L/2) = t k ipl um (L/2), the correctly 
scaled result is ifk (x) = tkf k mnl (x), where t k is yet un- 
known. We can readily express t k and from the bound- 
ary condition 



start with n = 1 



set the potential 
to be U new (n, x) 



set fc„ 



find numerically 
the smallest /c > fc m , n 
for which f n (k) = 



calculate for this k 
the wavefunction (z>(x) 



save k as fc n and 
ip(x) as (p n (x) 



stop if n = A/ ^ > 




No 



t fc (x) = e lkx + r fe e 



-ikx 



and from the flux continuity equation 



"2' 



(68) 



-f. (69) 



FIG. 19: Flowchart of a single Hartree-Fock iteration step as 
described in the text. 



The iteration step works as follows. We search for the 
Hartree-Fock energies e„ by solving the equation (fl9|) . 
Precisely, we solve the equation f n (k) — 0, where 



A similar procedure can be used to obtain the solution 
(f-k (x) together with the amplitudes t' k and r' k . 



f n {k) = Re{e- lkL /t k } - cos(27r0/^ o ) (71) 



Appendix B: Iteration steps of the Hartree-Fock 
calculation 

For numerical purposes it is useful-- to replace the 
potential r yS(x) + Uh(x) + Uf{u, x) in the Hartree-Fock 
equation (|35p by the potential 

U new (n,x) = f[-yS(x) + U H (x) + U F (n,x)} 

+ (1 - f)U old (n,x) . (70) 

In a given iteration step the equation ([35]) is solved with 
the potential U new (n, x), where U id is U new from the 
previous iteration step, and / < 1 is a properly chosen 
weight. In the first iteration step we set U new = ~/5(x). 



and tk is the transmission evaluated for the potential 
U n ew(n,x). Clearly, the equation f n (k) = is fulfilled 
for many different values of the variable k. Among these 
values we choose as physically correct only those k = 
k n which fulfill simultaneously two conditions. The first 
condition is the inequality e± < £2 < ■ • • < £n, where 
e%, £2, . • • £n are the energies of the N lowest levels, each 
of them given as e n = Y?k\/2m,. The second condition 
is that the wave functions <p n {%) of these N lowest levels 
are mutually orthogonal. The flowchart of the iteration 
step is shown in the figure 1191 

After finishing the iteration step we set the obtained 
<p n (z) hito the potentials Uh{x) and Up(n,x) and we 
obtain U new (n,x) for the following iteration step. We 
repeat the iteration steps until the energies e n do not 
change anymore. 
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We note that we calculate the transmission tk for the 
given U new (n,x) by means of the algorithm from Ap- 
pendix A. The wave function ip n (x) of the energy level 
£„ = h 2 k^/2m is calculated by using the procedure de- 
scribed in the last paragraph of Section ILA. 



Appendix C: Implementation of BBCI 

Here we outline the BBCI algorithm^. As mentioned 
in section II. E, we determine the single-particle basis 
ijjj(x) with the ladder of the single-particle-energy levels 
£j and we consider only the first N max levels by intro- 
ducing the upper energy cutoff as illustrated in the fig- 
ure [TJ Due to the cutoff, the expansion (|52")) contains the 
finite number of the Slater determinants Xn- This num- 
ber, equal to ( Nr ^ ans ^ by simple combinatorics, is usually 
huge. The question is how to asses the importance of all 
determinants and to omit those of them which 



(V) 



are unimportant. In principle, the importance criterion 
adopted in section II. E requires to compute the energy 



for all 



(Xn\H\x n/ — «" i N 
the obtained numerical values 



determinants and to order 
(Xn\H\Xn) increasingly. 



Since 



N 



is huge, such direct approach still con- 



sumes enormous amount of the computer memory and it 
is also time-consuming to order increasingly a series of 
^Nmax~j numbers. The bucket-brigade algorithm^ elim- 
inates these difficulties as follows. 

We use the single-particle states ipj (x) to generate the 
Slater determinants Xn via a recursion algorithm (Fig. 
|2T)| which is based on the sequence of buckets Sj^ with 
J = 0, • • • , N max and N = 0, • • • , N max , where N is the 
number of the particles in the bucket and J labels the re- 
cursion step J — ► J + 1. Each bucket Sj t N contains only 
the Slater determinants of N particles which occupy the 
single-particle states ipj(x) with j = 0, — , J — 1. In 
each recursion step J — > J + 1, we first expand the old 
buckets Sj^n by adding the Slater determinants from the 
buckets Sjjv— l) but with an extra particle created in 
the single-particle state j = J. After this expansion, 
the buckets Sj+i t N are truncated by selecting only the 
most important Slater determinants in accord with our 
importance criterion based on the chosen measure of im- 
portance, which is chosen as (Xn\H\x n ) in this paper. 
Owing to this algorithm our importance criterion is ap- 
plied only to the determinants in the bucket rather than 



to all 



(V) 



determinants. 



Appendix D: Matrix elements of electron-electron 
interaction 

If we expand the function f(x+y) into the Taylor series 
around x, we can write 



Buckets with Slater determinants X„ 

N=0 N=l N=2 JV=3 



step 
J 



copy 



m 




m m 



order by 
using 

and 
truncate 



step 
/+1 



\\\ 

W W \1 

> T T 



m 



WW] MMM 



m m 



BIB 



FIG. 20: Visualization of the recursion step J — > J + 1 dis- 
cussed in the text: Copy the remaining (after truncation) 
Slater determinants from the iV-particle bucket in step J into 
the iV-particle bucket in step J + 1. Copy the remaining (af- 
ter truncation) Slater determinants from the (N — 1) -particle 
bucket in step J into the iV-particle bucket in step J + 1, but, 
when copying, enlarge each Slater determinant by creating 
a particle in the state j = J. Truncate with the help of the 
measure (x„| H \xn) the Slater determinants of each bucket in 
the step J + 1. Go to the step J + 2. Stop the procedure after 
J reaches N max - The squares represent the occupied single- 
particle states j = 1, 2, . . . , J. The number of the particles in 
the bucket is labelled by N. 



L/2 

/oo 
dyV(y)f(x + y) = Y,Vnf {n) (x), (72) 



-L/2 



n=0 



where f^(x) = -gzxf{x) and 



L/2 

v n = ^ J dyV(y)y" 

-L/2 



(73) 



Setting the short range interaction V(y) = V$e \ y \/ d into 
([75]) we get 
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Vo 



L/2 



L/2 



J dye-W/V = ^[l + (-!)"] / d|/e-*/V 

o 

(74) 



-L/2 



and after integration per partes we obtain 



v n = v d[i + (-iy 



d n _ e -L/2d dm i L l 2 ) 



m=0 



(n — m) ! 



(75) 

Setting this v n back into (|72|) . reordering the summa- 
tions over m and n, identifying the particular Taylor se- 
ries of fW{x±L/2), and using the periodicity condition 
f( n >(x ± L) — p n '(x), we obtain the relation 



V a /3fS — 

-L/2 -L/2 

(79) 

into the form 



£ ^ /2tt<A' 

n=0 v 7 



^4; 



(80) 



where 



L/2 



dyV e-W*f( x + y ) 



-L/2 



2V dJ2 d2n \f (2n) {x)-e- L/2d f {2n) {x-L/2)\ . (76) 



n=0 



In the limit L 3> d this relation simplifies to 



L/2 

/ dyK )e -l^/ d /( a ; + 2/ )«2K ) dErf 2 "/ (2 " ) (^)- (77) 

J „ — n 



-L/2 



n=0 



Using the relation ([77)1 and substitution z = 2irx/L we 
can rewrite the matrix elements 



/3 7< 5 



(78) 



with aij a 7 s being one of the values { — 1, 0, 1} and with 



7T 

~. a x J 



(81) 

Here, the term with n = corresponds to the (5-function 
like electron-electron interaction. Due to the Pauli ex- 
clusion principle, this term does not contribute to the 
total energy and total current. The second term in 
the summation over n is proportional to d 3 . If we 
parametrize the exponential electron-electron interac- 
tion V(y) — Vbe~l y l/ d by the parameter ana an d ex- 
press Vo from that qrgi our matrix elements V%j ex- 
hibit two limiting cases. In the case Ak 2 F d 2 <§; 1, where 
Vq ~ ciRc/d 3 , the leading term in the summation over n 
is a d- independent constant. This is the reason why our 
many-body calculations give for small enough d the d- 
independent persistent current. In the case Ak 2 F d 2 ^> 1, 
where Vq ~ auc/d, the summation involves the terms 
proportional to d 2 , d 4 , etc, i.e., the matrix element is a 
complicated function of d. This is the reason why our cal- 
culations for large d give the persistent current depending 
on two parameters, d and Vo- 
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